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This work proposes a new approach to study transport properties of highly correlated local structures. The 
method, dubbed the Logarithmic Discretization Embedded Cluster Approximation (LDECA), consists of diago- 
nalizing a finite cluster containing the many-body terms of the Hamiltonian and embedding it into the rest of the 
system, combined with Wilson's idea of a logarithmic discretization of the representation of the Hamiltonian. 
The physics associated with both one embedded dot and a double-dot side-coupled to leads is discussed in detail. 
In the former case, the results perfectly agree with Bethe ansatz data, while in the latter, the physics obtained is 
framed in the conceptual background of a two-stage Kondo problem. A many-body formalism provides a solid 
theoretical foundation to the method. We argue that LDECA is well suited to study complicated problems such 
as transport through molecules or quantum dot structures with complex ground states. 

PACS numbers: 73.23.Hk, 72.15.Qm, 73.63.Kv 



I. INTRODUCTION 



The study of nanostructures has been motivated, on the one 
hand, by the potential appHcations in molecular electronics 
devices-i or in quantum computing^ and, on the other hand, by 
the search for a more profound understanding of fundamental 
many -body physics such as the Kondo effect. Experimentally, 
not only the existence of the Kondo effect in quantum dots^ 
or single-molecule transistors^ has been established, but it has 
also been demonstrated that nanostructures can be designed 
to produce more exotic phases such as multi-channel physics 
and thus, non-Fermi liquid behavior?^ On the theoretical side, 
while the single-impurity case is well understood by means 
of firmly established analytical^ and numerical methods, such 
as the Numerical Renormalization Group technique (NRG)<i 
the search for unconventional effects, non-equilibrium behav- 
ior, and the need to model complex real structures, such as 
molecules or multi-dot geometries, has triggered the develop- 
ment of alternative methodsi ^-^-'°-"''^-'^ For instance, the pro- 
cedure of exactly diagonalizing a finite cluster containing the 
many-body terms and embedding it into the rest of the system, 
the Embedded Cluster Approximation (ECA), has satisfacto- 
rily been used to study transport in nanoscopic structures in 
the last few years . '"'"I'^i'-^ Ideas similar to the embedded clus- 
ter method have been applied to the metal-insulator transition 
of the Hubbard model,iii^ 

Incorporating ideas from the Density Matrix Renormaliza- 
tion Group method (DMRG)'^ into NRG and vice versa has 
also resulted in substantial improvements in, e.g., the calcu- 
lation of dynamical properties'^ or time-evolution schemes,''' 
which now allows one to address problems previously out of 
reach for either method. In the same spirit, it is the objective 
of this paper to present the Logarithmic Discretization Em- 



bedded Cluster Approximation (LDECA) approach to study 
highly correlated electrons in nano-scale systems, combin- 
ing ECA with Wilson's idea of a logarithmic discretization of 
the conduction band.^ As one of our main results, we utilize 
many-body arguments to provide a solid theoretical justifica- 
tion of this formalism. Although the ECA method, due to the 
embedding process, is designed to analyze the infinite system, 
it produces results that depend on the cluster's size, which, in 
some cases, has led to controversial results.*2i2i LDECA not 
only successfully reduces these finite-size effects, but, more 
importantly, it also optimizes the description of the system in 
the vicinity of the Fermi level, allowing for the analysis of 
lower energy scales than accessible to ECA. 

To demonstrate the potential of the method, we focus on 
the physics of the Kondo effect in a single-dot, where we find 
excellent agreement with exact Bethe ansatz (BA) results. As 
there is a timely interest in more involved versions of Kondo 
physics, such as multi-channel situations,— SU(4),i22i21 as well 
as two-stage Kondo (TSK) effects , we further ap- 
ply LDECA to study a double-dot structure side-connected to 
leads. This system, with a subtle TSK ground state similar 
to the one studied in Ref. 2^ is an important testbed for our 
approach. Our results are encouraging, and we thus envision 
the successful future application of LDECA to more involved 
systems such as molecules adsorbed at metallic surfaces^*^ or 
dot structures with subtle ground states.— 

The plan of the paper is as follows. We first provide a dis- 
cussion of the theoretical foundation of the method in terms 
of diagrammatic perturbation theory in Sec. For the sake 
of a clear presentation, we choose to provide a pedagogical 
account of the theory, therefore the details will be given in 
an appendix (App. |A]i. Our results for the two systems, one 
embedded dot and a two dot model, are covered in Sees. |III] 
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binding lattice (thus, L = iV + 1, in the case of a single im- 
purity). In this cluster is where most of the many-body effects 
are expected to be confined. In a second step the cluster's 
Green function is computed with exact diagonalization, which 
then, in a last step, is embedded into the rest of the tight- 
binding lattice. ' 1- 12,13.20.23,3 1.32.33,34.35 jYie precise meaning 
of the embedding step is described below. 

The theoretical foundation of the method is outlined us- 
ing the Anderson single-impurity Hamiltonian describing a 
dot connected to a semi-infinite lead>^ The total Hamiltonian 
reads 



MB 



where 



and 



Hmb = U/2J2 



(1) 



(2) 



FIG. 1 : (Color online) (a) Schematic description of the steps needed 
in the implementation of the ECA and LDECA methods for the par- 
ticular case of one quantum dot coupled to a conduction band. The 
system is separated into two parts: the first one contains the many- 
body terms and the first few neighboring sites in the lead. This forms 
the cluster to be exactly diagonalized, with L = N + 1 sites. The 
second part is the rest of the lead, a semi-infinite tight-binding chain. 
A crucial step, which allows the simulation of Kondo physics, is the 
embedding of the cluster into the remaining part of the system. This 
step is performed through a Dyson equation, which amounts to a 
summation of an infinite family of Feynman diagrams arising from 
perturbation theory, (b) Diagrammatic expansion associated to the 
Dyson equation. The dressed propagators (in green, darkly shaded), 
which reestablish the artificially broken connection between sites N 
and R (through a hopping term V), are calculated as a function of 
the bare propagators (gray, lightly shaded). Note that the equation 
for the dressed propagator N+IJ does not have an independent bare 
term, since the bare propagator N+lj is zero. 



and Sec. lIVI respectively. For both models, we discuss the lo- 
cal density of states and the conductance as a function of gate 
potential. We close with a summary in Sec. |Vl 



II. THE "LOGARITHMIC DISCRETIZATION 
EMBEDDED CLUSTER APPROXIMATION" (LDECA) 

Similarly to the ECA method, LDECA is supposed to treat 
localized impurity systems that consist of a region with many- 
body interactions weakly coupled to non-interacting conduc- 
tion bands. The approach is based on the idea that the many- 
body effects of the impurities are local in character. With this 
in mind, we proceed in three steps: first, out of the complete 
system, one isolates a cluster with L sites that consists of the 
impurities plus their N nearest neighboring sites in the tight- 



Hi 



band 



^i(cLc,+i^ + 4_^i^Cj^) . (3) 



The first two terms of Ht represent the impurity, which has 
a diagonal energy, the gate potential Vg, and a Coulomb re- 
pulsion U in the Hamiltonian Hmb- The third term is the 
hybridization of the impurity with the band and finally, i^band 
represents the continuous spectrum, in this case modeled by 
a semi-infinite non-interacting chain, c]^ is a fermion cre- 
ation operator acting on site i, with a spin index a =ti- 
n-ia = is the particle density operator t' and are 

the hopping matrix elements between the dot and the leads 
and in the leads, respectively. A tight-binding band with a 
semi-elliptical density of states is obtained with the choice of 
U = 1. 

This problem can be treated within the framework of quan- 
tum perturbation theory. The standard many-body perturba- 
tion theory formulation generally considers the kinetic energy 
as the unperturbed Hamiltonian and the many-body terms as 
the perturbation. This permits the use of Wick's theorem to 
formulate a diagrammatic expansion for the propagators of 
the system. In our case, however, we adopt an opposite point 
of view. The unperturbed Hamiltonian consists of two parts, 
the isolated cluster, which includes the impurity and its neigh- 
borhood, and the rest of the system, as represented by the two 
dashed boxes in Fig. [T](a). The kinetic energy associated to 
the connection of these two subsystems is now considered to 
be the perturbation. This seems to be an appropriate starting 
point to describe a system where the many-body interactions 
are local, so that the cluster may contain most of the relevant 
physics we wish to describe. However, one faces several dif- 
ficulties in a theory of this kind. The most important one is 
the fact that, in this case. Wick's theorem is not valid and, as 
a consequence, it cannot be used to develop the diagrammatic 
expansion. However, perturbation theory provides us with a 
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way of proposing a locator-propagator diagrammatic expan- 
sion and establishing a criterion to sum up the most important 
families of diagrams (for details, see Appendix A). 

Therefore, following the strategy outlined above, the unper- 
turbed Hamiltonian Hq is given by: 
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where 



-^^cluster = ^a'^ "Ocr + Hhjb + (5) 
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(7) 



Figure [TJa) schematically represents the two parts of the 
system. Note that one of the internal connections of the lead, 
represented by a red line, labeled with ?lV — in the fig- 
ure, is artificially broken by this procedure and the two parts 
of the unperturbed Hamiltonian, i?ciuster ™d flrest> can be 
solved exactly. The ground state of i^duster "^^^^ ^ sites of 
the lead plus the impurity is obtained by using the Lanczos 
method;^ In addition, using a continued fraction scheme, the 
cluster Green functions at zero temperature are then evaluated. 
The Green functions for -ffrest are calculated exactly since it 
constitutes a one-body problem. 

To restore the artificially broken connection between sites 
A'' and + 1, the interaction between the cluster and the rest 
of the lead, 



h.C. 



(8) 



is taken as the perturbation in the many-body diagrammatic 
expansion for the Green functions. This step represents the 
embedding of the cluster into the rest of the system. 

For the sake of clarity, we restrict the discussion to the lo- 
cal diagonal Green function at the impurity site, while it is 
straightforward to calculate a non-diagonal Green function at 
two arbitrary sites i and j following the same prescriptions. 
To obtain the causal Green functions we follow the standard 
framework of an expansion in terms of Feynman diagrams?^ 
The causal Green function for the impurity site can be ob- 
tained from 



Goo^ait — t') 



T{c^„{t)cl{t')S{^)} 



(9) 



where, as usual, 5(oo) is the evolution operator and T is the 
time order operator. The mean values are calculated in the 
ground state of the unperturbed Hamiltonian Hq. 

The evolution operator S'(oo) is expanded in increasing or- 
ders of _ffp, which, when replaced in Eq. (|9]l, gives rise to a 
perturbation series for the Green function. 



The Green function of the system at the impurity, as dis- 
cussed in the appendix, can be written using a general Dyson 
equation, as schematically shown in Fig.[Tlb): 



•oo,.(a;) = G'il.iLo) + ^ G'^il{uj)Y.^{u)G,oA^) (10) 

i 



where i is restricted to be either or A^, w denotes frequency, 
and the self-energy Sf (ti^) is defined as 



(11) 



While Yj'^{uj) is a simple self-energy, represents an 

infinite expansion (see Eq. (|A19l l in the Appendix). It can 
only be calculated approximately, although this can be done 
in a systematic way by including terms in the expansion up 
to a certain order in U . Diagrams with a similar topological 
structure appear in the calculation of the one particle Green 
function for the Hubbard or Anderson impurity Hamiltonians 
treated in the thermodynamic limitfS In addition, in Ref. 
diagrammatic expansion for an interacting lattice in the strong 
coupling limit was used in order to include effects of long- 
range interactions beyond the exact diagonalization of a finite 
cluster 

The key approximation of LDECA is guided by a com- 
parison of the two contributions to the self-energy given in 
Eq. ( fTTI ). While Sq strongly depends on the size of the cluster 

through the non-diagonal Green function [Gq*)^ ^(w)]^, 
does not. This fact can be of great help in establishing a hi- 
erarchy between these two contributions to the self-energy. In 
order to achieve this, the applicability of this expansion is re- 
stricted to the vicinity of the Fermi energy, where we know the 
physics of the Kondo regime is contained. Following Wilson's 
logarithmic discretization of the lead's density of states,^ the 
Hamiltonian is rewritten by adopting hopping elements that 
depend on the site index i: 



(12) 



where A > 1 is a constant, i > 1 (i ~ being the posi- 
tion of the impurity), and we take i = 1 as the unit of en- 
ergy. Note that in the limit of A 1, the above expression 
for ti describes a semi-elliptical band, rather than the flat band 
commonly used in standard NRG calculations; however, close 
to the Fermi energy the two bands have the same low-energy 
physics. It is worth noting here that the above discussion ap- 
plies to both EGA and LDECA, with the exception that in 
EGA A is taken to be 1, implying that the band is not dis- 
cretized. 

The implications of this logarithmic discretization with 
respect to the contribution of the Hilbert space states are 
twofold: (i) near the dot, states of all energies are taken into 
account; (ii) of the states far from the dot, only those near 
the Fermi level are considered, while high energy ones are 
neglected.— Although by this procedure high energy scales 
are not well treated, it permits to accurately describe much 
smaller energy scales than it is possible with A = 1, for the 
same cluster size. This procedure is justified if the physics of 
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the problem depends only on states with energy close to the 
Fermi level, as it is the case in the Kondo effect, discussed in 
this paper 

As discussed in the appendix, around Eq. ( IA26b . from 
Eqs. dATg and (|AT9] i one realizes that T,% ~ A'^^"^) and 



f{N) A-(^-i'^/2^ such that 



fiN)X 



-{N-l){N/2-l) 



(13) 



The function f{N) is an intricate function of N which goes 
asymptotically to zero as N increases above the Kondo cloud 
length i^K, wich has a magnitude inversely proportional to the 
Kondo temperature^. This is the reason why even for A = 1 
(ECA), the self energy Sq will eventually become negligible 
in the limit of ^ ^i^- and can be disregarded for suffi- 
ciently large N. However, this is a length scale typically much 
larger than the value for which, according to Eq. ( fT3] l, the 
self energy Sq can be neglected in comparison to S^. Tak- 
ing A — \/2, for example, for a cluster with TV = 9, we get 
Sq/E^ ~ 10~^, reflecting the fact that, for a cluster size 
which permits diagonalization with a relatively modest nu- 
merical effort, the contribution to the self-energy can be re- 
duced to S^. This is a very favorable situation because Eq 
is a very complex object (see Appendix) that can be obtained 
only approximately, while E^ is very simple and can be cal- 
culated exactly. Therefore, here lies the key reason to intro- 
duce the logarithmic discretization into the procedure. 

Within this approximation, i.e., neglecting Eq, the embed- 
ding is carried out using Eq. ( IA18b and therefore, Eq. ( fTOt can 
be simplified to 



GooA^) = GlSA^)+G]^j^^V'gN+i,A^)GNo,A^) ■ (14) 

Note that the Green function of the semi-infinite linear 
chain at the site Nh-1, gN^i,a-{^), representing the leads, de- 
pends on the value of the parameter A. For A = 1, it is the 
Green function, g'"^, of a uniform semi-linear chain, given by 
g^'iu;) = {uj± ^^;2— 4i2)/(2t2). 

Before presenting results obtained with LDECA in the next 
two sections, we want to discuss some general aspects of the 
embedding procedure. If we were to study the low energy 
excitations by diagonalizing an undressed cluster without per- 
forming the embedding, the necessity of incorporating a large 
amount of states lying in the Kondo peak region would re- 
quire the diagonalization of a cluster of N sites such that 
< Tk, i.e., the energy scale associated with the bro- 
ken link V in Fig.[TJa) would have to be less than the Kondo 
temperature. In order to fulfill this condition, and at the same 
time choose a value of A that still adequately describes the 
neighborhood of the Fermi energy, the value of N would have 
to be such that the Lanczos diagonalization would become 
impractical due to the size of the Hilbert space. The embed- 
ding process solves this problem in a simple way by rendering 
the numerical diagonalization of a small cluster compatible 
with a correct description of the energy region immediately 
around the Fermi by allowing the contribution Eq((jj) to the 
self-energy to be disregarded. This approximation, even for 
A = 1, has shown to be surprisingly reliable to reproduce the 
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FIG. 2: (Color online) Sketch of the quantum dot geometries studied 
in this work, (a) Single quantum dot connected to two leads, (b) Two 
side-connected quantum dots studied in Sec. lIVI 



Kondo regime properties of various systems, showing suffi- 
ciently fast convergence with cluster size^^lls 



III. RESULTS: LDECA AND IQD 

In this section, LDECA is applied to study the conductance 
and the local density of states (LDOS) of a single quantum dot 
connected to two leads [see Fig|2ja)]. This case allows for a 
comparison of our results with an exact solution obtained from 
BA;^ Applying a standard basis transformation onto symmet- 
ric and anti-symmetric combinations of states (at sites located 
symmetrically with respect to the dot) the two-leads Hamilto- 
nian can be effectively written as having only one lead, ren- 
dering it identical to Eq. This shows that this example 
constitutes a one-channel Kondo problem4^ 



A. The local density of states 

We start with a discussion of the effect of the A- 
discretization on the LDOS. We expect that a larger density 
of poles close to the Fermi energy is induced by the dis- 
cretization, while fewer poles will be present away from Ep. 
The first aspect, the accumulation of poles close to E-p is ad- 
vantageous to properly describe Kondo physics. To still obtain 
a reasonable approximation to the LDOS away from the Ep, it 
turns out that it is preferable to use a a;-dependent broadening 
scheme, and we first detail this technical aspect. 

The dressed LDOS at the impurity is a collection of poles 
located at ojp, each one with its own weight Wp, given by the 
non-linear Dyson equation used in LDECA. As a consequence 
of the LDOS normalization, the weights satisfy J2p = 
2. In order to avoid distorting the LDOS curve through the 
artificially large separation of the poles away from the Fermi 
energy caused by the logarithmic discretization, and following 
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FIG. 3: (Color online) (a)-(c): Local density of states for A = 1, 2 
and 4, respectively, and a cluster size L = 8. Thin (red) lines rep- 
resent the LDOS as obtained directly from the LDECA procedure 
(ECA in the case of the top panel), while thick (black) lines are the 
LDOS for the same parameters, but after the broadening of the res- 
onances using (Eg. list has been applied, as explained in Sec lIII Al 
In the insets, we show a comparison of the LDOS of the bare clus- 
ter (thin (blue) line) and after embedding (thick (red) line) for small 
values of uj. All calculations done for (7 = 1.0 and U /ttF = 6.3. 



methods employed in NRGr^i^ we write the LDOS as a sum 
of logarithmic Gaussiansr^ 



exp 



(Inw-lnWp)' 



(15) 



where h is an arbitrary number that defines, together with Wp, 
the width of a pole located at cjp. We choose logarithmic 
Gaussians to represent the delta functions rather than the usual 
Gaussians or Lorentzians because this function is asymmetric 
with respect to ujy. This asymmetry, which effectively shifts 
the spectral weight of each pole to higher energies, compen- 
sates for the accumulation of poles at low energies in relation 
to higher energies, caused by the logarithmic discretization.'*'* 
As pointed out in Ref. [Sa, this procedure results in high- 
energy peaks that are slightly broader and asymmetric than 
in the case of the true LDOS. 

In Fig. [3] we show the LDOS at the impurity for the 
particle-hole symmetric situation = —U 12 and different 



values of A (thin red line) calculated with LDECA with a clus- 
ter of L = 8 sites. An imaginary part r\ = 0.001, common to 
all poles, was used for all curves that are obtained with a plain 
Lorentzian broadening of delta-functions (thin (red) lines in 
the main panels). The thick black lines show the LDOS using 
the logarithmic Gaussian broadening with h = 0.5. In this 
case, we obtain the characteristic LDOS for the Kondo prob- 
lem, consisting of a three-peak structure, with two of them 
located at = and u; = Vg + U and the third one, the 
Kondo resonance, located at the Fermi level Ep = 0. The 
notable difference between the LDOS for A = 1 and A > 1 is 
the sizeable naiTowing of the Kondo peak, in qualitative agree- 
ment with NRG^ (compare panels (a) with panels (b) and (c) 
in Fig.[3j. A more quantitative comparison with, e.g., NRG, 
will be presented elsewhere. We wish to draw the reader's at- 
tention to the inset of each panel, showing a comparison of 
the undressed LDOS [thin (blue) line] with the dressed LDOS 
[thick (red) line] . One can clearly see that the LDOS for the 
'bare' cluster (before embedding) vanishes at the Fermi en- 
ergy, while the LDOS after embedding is finite at = 0, 
corroborating the notion that the embedding step is crucial to 
capture Kondo physics. 

Figure 13 a) shows the LDOS for several values of the hy- 
bridization parameter F at the particle-hole symmetric point, 
and Fig. Ub) for a fixed ratio of U/ttT = 3.5 and several 
values of the gate potential. We define the hybridization as 
F = Tvpot''^, where po is the band density of states at the Fermi 
level. The top panel illustrates how the width of the Kondo 
resonance, i.e., Tk, decreases when F is reduced. In the bot- 
tom panel, we see how the LDOS for a fixed U/ttT = 3.5 
changes as the gate potential Vg is varied. Note that the Kondo 
resonance is pinned at Ep. In contrast, for |V^| > ?7 [dashed 
(green) curve, for Vg/U = —1-7] the quantum dot is dou- 
bly occupied and therefore there is no Kondo effect. In such a 
case, the LDOS has just one broad peak located at = Vg+U. 



B. The conductance 

Next, we demonstrate the effect of the A-discretization on 
the conductance. The conductance as a function of the gate 
potential Vg for a single embedded quantum dot is calculated 
by using the Keldysh formalism^^ In Fig. |5](a), we present 
the conductance for a fixed cluster size {L = 12) and different 
values of A, together with the exact value obtained by the BA^ 
For A = 1 (dot-dashed curve), there is a large discrepancy 
with the exact results: The conductance peak is too narrow, 
indicating that, in this case, the role played by the self-energy 
Sq cannot be neglected. However, for A = 2 (dotted curve), a 
value typically used in NRG calculations.fi our LDECA results 
substantially improve over the ECA ones [i.e., A = 1), and for 
values of A between 3 (dashed curve) and 4 (large-dots curve), 
the results accurately agree with BA. We have verified that 
LDECA reproduces BA results for U/T as large as tZ/F = 25. 

As discussed above, the effect of neglecting the self-energy 
Sq depends on both the value of A and the size of the cluster 
The dependence of the conductance on cluster size is shown 
in Fig.|5](b) for A = 4 and i = 4, 8, and 12. Beyond a cluster 
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FIG. 4: (Color online) (a) Local density of states for the particle- 
hole symmetric gate potential Vg = — [7/2 for U = 1.0, A = 2, and 
different values of U/ttF, going from the Kondo regime (U/tvF = 
6.3) to the intermediate valence regime {U/ttF — 1.5). (b) LDOS 
for three different values of the gate potential for U/ttT = 3.5 and 
A = 2. The inset shows a comparison of results for Vg/U — —0.24 
between A = 2 [solid (black) curve, same as in the main panel] and 
A = 1 [dashed (blue) curve]. This illustrates the better 'pinning' of 
the Kondo peak to the Fermi energy, achieved with A > 1. 




FIG. 5: (Color online) Comparison of LDECA and BA. Conductance 
vs. Vg for U = 0.96 and t' = 0.2. (a) LDECA for L = 12 and the 
values of A indicated, (b) LDECA for L = 4, 8, and 12, with A = 4. 
The BA result is the solid black line in both (a) and (b).= 



Vg /U = —0.24 curve presented in panel (b), which was calcu- 
lated with A = 2. The dashed (blue) curve has been calculated 
with the same parameters, but for A = 1. The comparison 
clearly shows that the A = 2 result has more spectral weight 
at the Fermi energy than the A = 1 result. This increase of the 
spectral weight in the LDOS is at the heart of the improvement 
achieved for the conductance by using the band discretization. 



IV. RESULTS : LDECA AND A TWO-STAGE KONDO 
SYSTEM 



size Lc{X), the conductance is almost independent of L. For 
instance, for A = 4, results for i > 8 are indistinguishable. 
This characteristic length Lc{X) decreases as A increases. As 
A controls the extension of the neighborhood of the Fermi en- 
ergy that is accurately described, i.e., the larger the value of 
A, the smaller this region is, a compromise has to be found 
between the size of the cluster and the extension of the energy 
region around the Fermi energy that needs to be accurately de- 
scribed. Obviously, this depends on the model and the prop- 
erty being analyzed. The important point to be emphasized 
is: the results in Fig. |5]show that, with a value of A similar to 
the one widely used by the NRG community, it is possible to 
reproduce the exact results using a cluster size accessible to 
the Lanczos algorithm. 

It is also interesting to note that the improvement of the 
conductance results for A > 1 as compared to A = 1 are 
associated with a better 'pinning' of the Kondo peak to the 
Fermi energy. This can be partially inferred from the LDOS 
results shown in Fig. |4jb), where the two solid curves have 
the Kondo peak pinned at the Fermi energy. This statement 
can be made more quantitative by considering the results in 
the inset of Fig. SJb), showing a comparison between A = 1 
and 2. In that inset, the solid (black) curve is an enlarged 
view of the LDOS at the vicinity of the Fermi energy for the 



A. Overview: Regimes of the model 

We next analyze a system composed of a double-dot side- 
connected to a lead. The inter-dot and dot-lead connections 
are given by the matrix elements t" and t', respectively, as 
sketched in Fig. |2b). The transformation to symmetric and 
antisymmetric states is applied, since we again deal with a 
one-channel Kondo problem. After performing that transfor- 
mation, the Hamiltonian is given by: 

= Vg^Uda + U/2^ndanda 

+t'\/2^(4^^cia + cl^Cd^a) + i?band(16) 

OO 

i—ltT 

where we use the ti as given in Eq. ( fT2] i. and d = di,d2, 
labeling dot 1 and dot 2, respectively. 

The transport properties of this two-dot system can be ex- 
pected to be controlled by the interplay between the Kondo 
effect and the antiferromagnetic inter-dot correlation, and by 
the interference arising from the two distinct paths available 
to the electrons: visiting or bypassing the dotsi^ 
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As previously discussed in the literature, systems similar to 
the one depicted in Fig. Hfb), such as, for instance, the so- 
called T-configuration,i2^ exhibit two distinct regimes de- 
pending on the ratio t" jt': (i) when t" 3> t' , one is in the 
molecular regime and (ii) for t" ^ t' , the system crosses over 
into the TSK regime. It is important to realize that indepen- 
dently of t" , we expect perfect conductance at = —Ujl. 
Indeed, at the particle-hole symmetric point the dots always 
form a singlet, which is of different nature though, depending 
on the ratio t" jt' , as explained below. 

In the molecular regime, on the one hand, both dots act as 
a single entity, in a way that, as a function of the gate poten- 
tial, whenever an overall finite magnetic moment is located in 
the structure, the system exhibits a single-stage Kondo effect. 
In this regime of t" ^ the system essentially behaves as a 
single-dot with the two relevant levels separated by a large en- 
ergy. On the other hand, in the limit of a small t" jt' , such that 
the effective antiferromagnetic spin-spin interaction between 
the dots satisfies J=A{t"Y /U < Tk, the system is in the two- 
stage Kondo regime, which is characterized by a new energy 
scale To related to dot 1, much lower than the Kondo temper- 
ature Tk associated to dot 2. For T < To, very low energy 
physics is involved^^^i^ difficult to be captured by numerical 
methods such as standard ECA or DMRG. Yet, as shown be- 
low, a correct result for the conductance can be obtained from 
LDECA. 

Note that TSK behavior may manifest itself both as a func- 
tion of temperature at a fixed gate potential2§^ and as a function 
of gate potential at a fixed temperature;^ As our method is a 
zero-temperature one, we will here focus on the gate potential 
dependence of the conductance and other quantities. We will 
argue that as one starts from the empty orbital regime, first a 
single Kondo effect emerges at > 0, which causes a sup- 
pression of the conductance. As the gate potential is further 
tuned towards Vg = —U/2, the magnetic moment of dot I is 
eventually Kondo-screened as well through the quasi-particles 
of the composite system of dot 2 and the lead. This gives rise 
to the aforementioned singlet between dots 1 and 2, which 
leads to perfect conductance atVg = ~U/2. 

The plan of this section Sec.|lV]is thus the following. We 
first demonstrate in Sec. IIVB] that indeed, a singlet is formed 
between the two dots, independently of t" /t'. These results 
are obtained with both DMRG and a diagonalization of the 
bare clusters, before the embedding process. Second, we com- 
pute the conductance and charge as a function of gate potential 
and discuss these results both in the molecular and the TSK 
regimes in Sees. IIV Cl and UVDl respectively. We further aim 
at illustrating how the properties of the system change as it 
crosses over from the molecular regime into the TSK regime. 
Finally, at particle hole symmetry {Vg = — C//2), we present 
LDECA results for the LDOS at the dots and the conductance 
as a function of t" /t' . As a key result, we demonstrate that 
using the discretization of the band, LDECA produces perfect 
conductance down to very low values of t" /t', a result which 
was previously out of reach for ECA (A — 1). 
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FIG. 6: (Color online) Spin-spin correlations calculated by Lanczos 
and DMRG as a function of Vg for U = 1.0 and t' = 0.3. (a) t" = 
2.0 (molecular regime) and (b) t" — 0.05 (TSK regime). Lanczos 
with L = 11 and A = 6, (Si ■ S2) (solid) and (S2 ■ Sc) (dashed) 
black lines; DMRG with 196 sites at half filling of the full system 
and A = 1, (Si ■ S2) (circles) and (S2 • Sc) (squares). 



B. Spin-spin correlations 

For the present model, we now establish the presence of a 
strong antiferromagnetic correlation between the dots by ana- 
lyzing the spin-spin correlations as a function of V^, presented 
in Fig. |6]for both t" = 2.0 (upper panel) and t" = 0.05 (lower 
panel), with t' = 0.3 in both cases. Results are for large clus- 
ters with 196 sites and A = 1, obtained with DMRG, and 
also for L = 9 and A = 2 (upper panel), and for L = 11 
and A = 6 (lower panel), using a Lanczos diagonalization 
procedure. Both DMRG and Lanczos calculations were done 
without embedding. 

In the case of the molecular regime [t" = 2, Fig.|6|a)], the 
inter-dot correlation (Si -82) is large for —2.5 < Vg/U < 1.5. 
Therefore, we expect a perfect conductance in that window, 
and a Kondo anti-resonance to appear at Vg/U > 1.5 (see 
Fig.lHa), below). For the smaller t" = 0.05 [Fig. Iltb)], the 
antiferromagnetic correlation between the dots is dominant in 
the window —U<Vg<Q, indicating the formation of a sin- 
glet. While the inter-dot spin correlation (Si • S2) is large at 
the electron-hole symmetric point Vg = —U/2, the antifer- 
romagnetic correlation (S2 • S,;), although small, is not zero. 
Note that site c is adjacent to dot 2, see Fig. [2b). For instance. 



8 




-0.1 
-0.2 

A -03 
in 

t/T -0.4 
V 

-0.5 
-0.6 
-0.7 





<S,S,> i.= 6, 11 sites 

B--n <SjS^>DMRG 196 sites .qqi 

<S,S,> ^ = 6, 1 1 sites 

o- -o <S,S,> DMRG 196 sites '^^-O 02 

' - m 

V 

-0.03 
-0.04, 









0.5 



t'Vt' 



1.5 



FIG. 7: (Color online) Spin-spin correlations as a function of t" jt' 
for U = 1.0, t' = 0.3 and Vg = -f//2. See legend and Fig.|6]for 
the description of the lines. 



(S2 • Sc) takes a maximum at w in the case of t" = 0.05, 
indicative of the single-stage Kondo effect that is observed in 
that gate potential region [see Fig. fTTT a)]. 

We now study both correlations as a function of t" /t' at 
Vg = —U/2, which is displayed in Fig. Q An important 
observation is that (S2 • Sc) increases in magnitude as t" is 
reduced, as is shown in the inset of Fig. |7] This fact indi- 
cates the subtle existence of a Kondo-like ground state, which 
is strengthened when t" is reduced. In addition, the antifer- 
romagnetic inter-dot correlation, presented in the main panel 
of Fig. |7] also increases when t" is reduced, taking values 
as large as (Si • S2) ~ -0.7, for t"/t' > 1/10. This is 
surprising, since the interplay of these two correlations has 
a behavior opposite to other known systems, such as heavy 
fermions near a quantum phase transition^ or embedded two- 
dot nanostructures.— It reflects the existence of an inter-dot 
singlet for all values of t">Q. However, the nature of the sin- 
glet for small t" /t' is different from that for the large t" /t' 
regime. While in the latter case the singlet, which is caused 
by the direct interaction between the dots, destroys the Kondo 
regime, in the former case it is enhanced by the Kondo spin 
correlation with the intervention of the conduction electrons, 
as shown by the fact that the inter-dot and the Kondo spin 
correlations increase when t" is reduced. These results illus- 
trate the characteristics of a TSK state. The comparison of 
data from large clusters (196 sites, A = 1) and short ones 
(L = 11, A = 6) yields a convincing agreement for the spin 
correlations, especially for (Si • S2), governing the inter-dot 
singlet. This agreement indicates that the spin-spin correla- 
tions are, so to speak, localized objects. The embedding pro- 
cess is not as important to calculate static properties as it is 
for the conductance, which we shall see later. Still, the van- 
ishing of (Si • S2) at very small t" /t' on the smaller cluster 
reflects that, in this particular limit, the embedding is crucial 
to overcome this finite-size effect. An important point that we 
want to emphasize here is that the molecular and TSK results 
suggest perfect conductance at = —U/2 for any nonzero 
t". 
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FIG. 8: (Color online) Schematic representation of the molecular 
regime, (a) When t" 3> t' , the two dots form a single entity, rep- 
resented by the shaded box. (b) - (d) Different regimes at different 
values of gate potential: (b) Kondo effect for the bonding orbital 
(zero conductance, because of the back-scattering density of states 
at the Fermi energy Sf), (c) Perfect conductance at Vg — —(7/2 
(no back-scattering density of states at Ef), (d) Kondo effect for the 
anti-bonding orbital (zero conductance). 



C. Molecular regime 

We proceed with an analysis of the conductance for the so- 
called 'molecular' regime. In Fig. [HJa), we schematically il- 
lustrate what happens for t" /t' ^ 1.0, i.e., when the inde- 
pendent dots are 'locked' into 'molecular' bonding and anti- 
bonding orbitals, separated by a large energy, proportional to 
t". The two dots now behave as a single structure, represented 
by the dashed square box, side-connected to the leads. The ef- 
fect of these molecular orbitals over the conductance through 
the leads, as the gate potential Vg varies, is pictured in panels 
(b) to (d), where now the bonding and anti-bonding orbitals 
are depicted inside a quantum well. In panel (b), the bond- 
ing orbital is in the Kondo regime. As the double-dot struc- 
ture is side-connected to the leads, the conduction electrons 
are back-scattered, resulting in zero conductance^^ Panel (c) 
shows that, at the particle-hole symmetric point Vg = —U/2, 
the bonding orbital is doubly occupied, lying below the Fermi 
energy Ep, and the anti-bonding orbital, lying above Ep, is 
empty. Therefore, in this case, the double-dot structure creates 
no back-scattering density of states at Ep [as schematically 
indicated in the panel (c)], resulting in perfect conductance. 
Finally, panel (d) displays the corresponding Kondo effect for 
the anti-bonding orbital, also resulting in zero conductance. 

Figure |9] shows the conductance and the charge vs. gate 
voltage in the molecular regime as obtained with LDECA for 
t" = 2.0, t' = 0.3, and [/ = 1 on a cluster with L = 11 sites 
and A = 2. In this plot, one observes the two Fano-Kondo 
anti -resonances, with an approximate width of U/2, originat- 
ing from two 'molecular' levels separated by w 2 x t" + U/2. 
It is important to emphasize that in the molecular regime the 
two dots behave as a unique entity, providing an extra lateral 
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FIG. 9: (Color online) (a) Conductance and charge as a function of 
Vf, for the molecular regime (U = 1.0, t' = 0.3, and t" = 2.0) 
for the two side-connected quantum dots. Notice that each feature of 
the conductance is associated to the respective diagram introduced 
in Fig. [8] Panel (b) displays the LDOS for the bonding state (B) 
for V^/U = 1.78 (the anti-bonding state (AB) has an identical peak 
for Vg/C/ = -2.78). Panel (c) displays the LDOS for the bonding 
(solid (black) line) and anti-bonding (dashed (red) line) states at the 
particle-hole symmetric point Vg — —U/2. A = 2 for all calcula- 
tions. 



FIG. 10: (Color online) Schematic representation of the TSK regime, 
obtained when t" <^ t' . (a) Sketch of the first Kondo stage, when 
the two dots are occupied by a single electron, which forms a Kondo 
singlet with the conduction electrons, represented by the vertical an- 
tiparallel arrows (one in the dots, the other in the band). The Kondo 
coupling, which mediates the formation of the many-body singlet, is 
represented by the (yellow) lightly shaded oval. The density of states 
created by the Kondo effect in dot 2 at the Fermi energy induces 
backward scattering of the conduction electrons going through the 
leads (horizontal antiparallel arrows). This results in zero conduc- 
tance, (b) AlVg = —U /2, when the dots have total charge n = 2, 
the spins in the dots are locked into a singlet (darker shaded oval). 
This singlet is mediated through two Kondo effects, represented by 
the underlying lighter (yellow) oval. The combination of these two 
Kondo effects, in contrast to the case in panel (a), removes density of 
states from dot 2 exactly at the Fermi energy, suppressing the back- 
scattering and resulting in perfect conductance. This is represented 
by the rightward arrow 'piercing' the Kondo coupling oval. 



path for the electrons to traverse when visiting the Kondo peak 
derived from the molecular orbital. This gives rise to the Fano 
antiresonance in the conductance appearing in Fig. |9ja). As 
expected, both dots are charged almost simultaneously, as one 
can see in Fig. |9ja), with a dashed line for dot 1 and a solid 
thin line for dot 2. In the top of Fig. |9] the corresponding 
potential wells described in Fig.[8]are displayed. 

To illustrate the idea of a 'molecular orbital Kondo effect', 
in the lower left panel, we display the LDOS associated with 
the molecular bonding orbital formed with the two dots. This 
LDOS is calculated at the positive gate potential at which the 
conductance is zero, which turns out to be at Vg « t" — U / 4, 
as expected. We find the 'molecular' Kondo peak at the Fermi 
energy, as well as the broadened e and e + U levels, where the 
renormalized intra-orbital Coulomb repulsion U = U /2 can 
be obtained by rewriting the dot Coulomb repulsion in the ba- 
sis of the bonding and anti-bonding orbitals. A similar result 
(not shown) is found for the LDOS of the anf/bonding orbital 
at the Vg value for which the second Kondo-Fano resonance 
occurs. The LDOS for each quantum dot (not shown) also ex- 
hibits a Kondo peak. Indeed, since the two dots equally partic- 
ipate in the molecular Kondo effect, their LDOS are quite sim- 
ilar to each other, and qualitatively similar to what is shown 



in Fig.|9jb). Finally, in Fig.|9|c), we show the LDOS for both 
the bonding and anti-bonding orbitals for Vg = —U/2, where 
clearly the Kondo peak is absent and the two orbital levels are 
separated by about 2t" . 

D. Two-stage Kondo regime 

Figure[TO]schematically depicts a much more subtle regime 
than that of Fig. [8j the TSK regime. One enters into this 
regime when t" <^ t', where now the connection of dot 2 
with the leads is much stronger than the inter-dot connection. 
Here, the concept of bonding and anti-bonding orbitals does 
not apply, since each dot feels the interaction with the con- 
duction electrons differently, the crucial point being that dot 1 
interacts with the Fermi sea through dot 2. 

In this and the next section, we will present evidence that 
our LDECA results are perfectly consistent with the notion of 
TSK behavior. As a guidance to interpreting the numerical 
results, this behavior can be schematically described as fol- 
lows. In Fig. [TOta). where the gate potential Vg is such that 
the charge occupancy of the two dots is 1, i.e., ni + n2 — 1, 
a Kondo effect develops, represented by the oval shape with 
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FIG. 1 1 : (Color online) (a) Conductance and (b) charge in the dots 
as a function of Vg for values of t" < t' (U = 1.0, t' — 0.3, and 
values of t" as indicated). See text for discussion. All calculations 
are for A — 2.0, except for t" = 0.125, where A = y/6. Notice 
that the conductance and charge features at Vg « and Vg = —U/2 
were associated to the schematic diagrams introduced in Fig. 1101 



dashed borders, resulting in back-scattering and a vanish- 
ing conductance, as indicated by the horizontal arrows. The 
Kondo effect involves a magnetic moment located in the dots, 
which is screened by the conduction electrons, indicated by 
the arrow on the dot and an antiparallel one in the band. 
In Fig. [Tolb). depicting the situation for Vg = —U/2, the 
two dots are each singly occupied and a strong singlet forms 
between them, represented by the darker oval shape with a 
solid border. Although in this regime the LDOS of dot 2 is 
exactly zero at the Fermi energy, therefore suppressing the 
back-scattering and restoring perfect conductance (this is rep- 
resented by the arrow 'piercing' the lighter shaded oval), it 
does not eliminate the Kondo spin-spin correlation between 
the spin of dot 2 and the conduction spins. On the one hand, 
the conduction electrons do not see the two dot system as a 
unique entity: Indeed, they recognize dot 2 as a separate ob- 
ject to which their spins correlate. On the other hand, the spin 
of dot 1 sees the rest of the system as a whole, and Kondo cor- 
relates with the spin of dot 2, thus creating the two-dot singlet 
state. In reality, the singlet is a many-body effect involving 
the conduction electrons and it is composed of two consecu- 
tive Kondo effects (represented here by the underlying lightly 
shaded oval). 

In Fig. [TTl we show the LDECA conductance [panel (a)] 
and the charge in each dot [panel (b)] as a function of Vg 
for much lower values of t" than in the previous molecular 
regime, namely t" = 0.2, 0.15, and 0.125. Let us first discuss 
the charge, as an example of a quantity that exhibits a qualita- 
tively different behavior for high and low values of t" / 1' , with 
the two dots behaving more independently as t" 0. 

In the small t" /t' regime, dot 2 is charged first upon ap- 
proaching = 0, and only when it has a substantial amount 
of charge dot 1 starts to be charged as well [see Fig. fTTT b)]. In 
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FIG. 12: (Color online) Conductance vs. t" /t' for U = 1.0, t' = 
0.3, Vg = —U /2. Results are for L — 11 and A = 1 (black circles), 
A = 2 (red squares), A = 4 (green diamonds) and A = 6 (blue 
triangles). In the inset, we show t"(A)/t' for L = 11 (black circles) 
and L — 5 (red squares). 



addition, around Vg ~ where the minimum in the conduc- 
tance occurs due to the single-stage Kondo effect, the curve 
for the charge of dot 2 features a much more well defined 
plateau than that for dot 1. This indicates that, in contrast 
to the molecular regime [see charge behavior in Fig. |9ja)], 
the two dots now start to have a qualitatively different par- 
ticipation in the Kondo effect, suggesting that at this much 
lower value of t" /t' and at « 0, one starts to see the 
emergence of the first stage of the TSK regime. (Notice that 
DMRG results for the charge and the total spin as a function 
of Vg for large clusters (not shown) agree with the LDECA 
picture just described). In addition, the width of the Kondo 
anti-resonance seen in Fig. [TTl a) is now substantially smaller 
than U/2, which is the typical value found in the molecular 
regime, see Fig. |9ja), although much larger than the intrin- 
sic width of the dots' resonance states. Finally, as one ap- 
proaches Vg = —U/2, and each dot now has one electron, the 
second stage of the TSK is reached. In this regime, through 
the mediation of the conduction electrons, the interdot-singlet 
is formed and, although it shows a two-Kondo peak structure, 
the LDOS of dot 2 is zero at the Fermi energy [see Fig. [T3T b)l. 
As a consequence, the system exhibits perfect conductance 
[see Fig. fTTT a)]. 



E. Conductance and LDOS at Vg = -U/2 

Next, we discuss the conductance at the particIe-hoIe sym- 
metric point as a function of t" to show that LDECA correctly 
captures the low-energy physics down to small values of t". 

Figure[T2]displays the conductance vs. t" for various values 
of A and L ~ 11. We recall that, as exemplified in Figs.|9]and 
[TTl the conductance at = —U/2 for any t" > should be 
G = Go- We study the electron-hole symmetry situation, as 
it is the most difficult point to be correctly described, having 
the lowest Kondo temperature for the set of parameters taken. 
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The suppression of the conductance for small values of t" /t' 
shown in Fig.[T2]is caused by finite-size effects which obscure 
the second stage Kondo effect. In this specific case, we find 
the tendency of a strong suppression of spin fluctuations in 
dot 1 as the system approaches half-filling, {Vg = —U/2). 
This behavior at A = 1 is similar to other models discussed 
in detail in Ref. [HI Figure [12] suggests that the finite-size 
dependence of the conductance for A = 1 (circles) is quite 
severe, as it starts to manifest itself ait" /t' k, 0.5. However, 
it is also evident that by increasing A the situation improves 
markedly, which is the main message to be taken from this 
figure. 

From the curves for each different A we can extract a char- 
acteristic inter-dot coupling t"{\), satisfying G{t")/GQ = 
1/2, below which the conductance rapidly approaches zero. 
The dependence of t" on A~^ for two different values of L is 
shown in the inset. t"{X) tends to zero for values of A that de- 
crease with increasing cluster size. Therefore, at Vg=—U/2, 
G/Go^l wheni" ^ 0. 

To further demonstrate the difference between the molec- 
ular and the TSK regimes at = ~U/2, Fig. [T3] shows a 
comparison between the LDOS for two widely different val- 
ues of t"/t'. In Fig. [Bja), we show the LDOS at dot 2 for 
t" = 2.0, t' = 0.3, U = 1.0, and A = 2. These are the same 
parameters as the ones used in Fig. |9|c), where the bonding 
and anti -bonding orbitals where shown. Figure [T3l a) unveils 
why the conductance in the molecular regime is G = Go at 
Vg = —U/2 [see Fig. |9ja)], as dot 2 has a vanishing den- 
sity of states in a wide energy region around the Fermi level. 
Once there is no back-scattering density of states at the Fermi 
energy, the conductance is perfect. 

On the other hand, in Fig. [Olb), LDOS results for dot 2 
are depicted for t" = 5^^, and A — V6, and the same val- 
ues of t' and U as in Fig. [TSl a). Again, the density of states 
at the Fermi energy vanishes. However, in this case, close to 
the Fermi energy, we find two sharp features, suggestive of 
a Kondo peak split in two. To substantiate this picture, the 
dashed (red) curve in Fig.[T3jb) shows the Kondo peak that is 
present when t" = 0, i.e., when dot 1 is effectively removed. 
It is the presence of dot 1, interacting with the rest of the sys- 
tem through dot 2, that gives rise to the TSK regime, reflected 
in the LDOS of dot 2 as an antiresonance in the middle of its 
Kondo peak. 

In summary, the fact that the Kondo regime of dot 1 is medi- 
ated by the Kondo state of dot 2 explains the surprising result 
of a perfect conductance at Vg = —U/2 in the TSK regime. 
This mechanism reduces the LDOS at the Fermi level of dot 2 
to zero as shown in Fig. [T3l b). ehminating an alternative path 
for the circulating electrons and hence any destructive inter- 
ferences. In this regime, the electrons at the dots form a spin 
singlet, even at small t" /t'. This subtle effect and its con- 
sequences on the conductance are well captured by LDECA. 
It is then clearly shown in Figs. [T2l and [T3] that the logarith- 
mic discretization of the band, combined with the embedding 
process, provides reliable results in a wide parameter range. 
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FIG. 13: (Color online) (a) LDOS at dot 2 for i = 11, t" = 2, 
t' = 0.3, U — 1, Vg = —U/2, and A = 2 (same parameters as 
the ones in Fig.[9jc), but here we show the LDOS of dot 2 only). 
Notice the large region of vanishing LDOS around the Fermi energy 
{uj — 0). (b) Same as in (a), but now for t" — 5x 10~^ and A — 
(solid (black) line). Notice that, as in (a), the LDOS vanishes at the 
Fermi energy, but now only in a very narrow interval. A comparison 
of the width of the double-peak structure at low cu with the LDOS 
results for t" = (dashed (red) line), i.e., for a single side-connected 
dot, suggests that the double-peak structure is a split Kondo peak. An 
enlarged view of the double-peak structure is seen in the inset. 



V. SUMMARY 

In this paper, we developed a formalism to study local and 
highly coiTelated electrons that combines the numerical sim- 
plicity of the ECA method with Wilson's idea of a logarith- 
mic discretization of the non-interacting band. A diagram- 
matic expansion that provides a solid theoretical basis for 
the method was also discussed. Applied to a one-impurity 
problem, LDECA yields an excellent agreement with BA re- 
sults. In addition, following the same procedure used in NRG 
to broaden the LDOS, a perfect agreement was found with 
the accepted results for the LDOS of the Anderson impu- 
rity model. In the case of a double-dot side-connected to a 
lead, and at small t", the contrast between the A=l (ECA) and 
A > 1 (LDECA) results exemplifies the power of the loga- 
rithmic discretization: the low-energy physics associated to 
the TSK regime is correctly unveiled, as LDECA provides an 
accurate description of the physics close to the Fermi energy. 

The main advantage of LDECA is its great flexibility, which 
allows the incorporation of other degrees of freedom, such 
as localized phonons or photons.'*^ The restrictions imposed 
by the Lanczos method can be overcome by using DMRGr^ 
allowing one to use larger systems and to study more involved 
problems. Related efforts are in progress. 

We conclude that the LDECA method can be applied 
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to complex problems, including molecules adsorbed at a 
metallic surface22 and sophisticated topologies of quantum 
dots, displaying exotic Kondo regimes, such as, for example, 
non-Fermi liquid behavior, two-channel Kondo effect and the 
physics associated to SU(N) systems. 
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It is important to emphasize, as shown in Eq. ( IA3I ). that the 
conservation of charge of the unperturbed subsystem restricts 
the expansion to even orders in V. The undressed Green func- 
tion appearing in the equation is defined as 

50,7V,JV+l,cr,o-i,o-2 (^1 ^^27 i') — 

('^{co<T(i)c]v^i(il)cAr+l<Ti(il)c]v+i^,(i2)cAra2(t2)cS^(t')})^ 

(A4) 

with an obvious generalization for the undressed Green func- 
tion of other orders. Calculating terms of all orders in V in the 
expansion, Eq. (IA3l l, the local Green function can be written 
as 

Goo,.(t-i') = 41(i-^')+G^ol.(^-i')+4ol(i-^')+- • • ■ 

(A5) 

As the operators belonging to the two different unperturbed 
parts of the system are, in this ground state, decoupled from 
each other since there is no connection between the cluster 



APPENDIX A: DIAGRAMMATIC EXPANSION 

In this appendix we develop the diagrammatic expansion of 
the one particle Green function at the impurity site given by 

(T{coAt)cl{t')S{^)}) 
(5((X)))o 

where, as usual, 5(00) is the evolution operator and T is the 
time order operator. The mean values are calculated in the 
ground state of the unperturbed Hamiltonian Hq, given by 
Eq. (HJi, restricting our discussion to zero temperature. 

The evolution operator 5* (00) is expanded in increasing 
orders of the perturbing term Hp, which, when inserted in 
Eq. ( lAlb . gives rise to a perturbation series for the Green func- 
tion. It can be written as 



(A2) 



(A3) 



and the rest of the leads, Eq. ( IA4I ) results in 

gO,N,N+l,(T,a-i,a-2 (^j ^1,^2,^') 

= 90,N,a,criit,ti,t2,t')gN^i^criitl,t2)Saia2 (A6) 

where the spin conservation imposes the condition cti = (T2 
and 

gO,N,a,aiit, ti,t2,t') = (T {co^^{t)c'^ ^^^{ti)cN,aiit2)4)^^{t')^^ 

(A7a) 

gN+l,ai{tl,t2) = (T'{cN + l.aiitl)c'N+l^^^{t2)}'j^ •(A7b) 

According to Eq. iA6i . the expansion Eq. iA3\ is formally a 
locator-propagator expansion^ where the locators correspond 
to the unperturbed sub-systems Green functions and the prop- 
agator turns out to be the one connecting them. 

The Green function go^N.^.a-i can be diagrammatically rep- 



^(°°) = E ^T^r^ / ■■■/ (T{co.WiJp(tl)...i?p(t„)}ct^(i')) dt...dtr. 

,7^1 n\{S{oG)) Jto Jto^ /o 
Substituting this expression into Eq. (lAlb . the local Green function is given by 

Goo,a{t ^ t') = g0.a{t — t') + J go^N,N+l,a,cri,cr2{t7tlTt2,t')dtidt2 

(71(72 

+V'^ ^ g0,N,N+l,<T,cri,iT2,cr3,cr4itjtl7t2,t3,t4,t')dtldt2dt3dt4 + 
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resented by 

iJi \ai 

t a 



l-O 

t' 



r 



\coa{t)c^Na^ (^l)cWai (^2)4^ (^') } 



and 



gN+l.aitlih) 



ti a to 



The zeroth-order contribution to the Green function, i.e., the 
solution of the problem for \^ = 0, is represented in terms of 
diagrams as 



= O-l l-O 

t <7 t' 

(A8) 

The Green function is exactly obtained by calculating the 
ground state of the cluster using the Lanczos method. Al- 
though it includes the many-body interaction and its effects 
within the cluster, it is the undressed Green function with re- 
spect to the expansion given by Eq. ( IA5b . 

To second order in perturbation theory, the contribution to 



'oo.a(*-^')is 



Gg>it-t') 



(5(oo))o 

jgaMMai (i, ^1,^2, t')9N+l,ai {ti , t2)dtidt2 (A9) 



and can be diagrammatically represented as 



CTl 

O-C 

t <J 



t-2 
=1-0 



(AlO) 



where an integration over each internal time t\ and is im- 
plied and the sum over a\ needs to be taken. 

Regarding the calculation of the many-particle Green func- 
tions at the lead sites outside the cluster. Wick's theorem can 
be used since this part of the system is represented by a one- 
body Hamiltonian. In this case, it is clear that for the fourth 
order we obtain the diagram 



1 



G''^lAt-i!) = ^ 
t 



0-1 



t2 h 

(Tl (72 



0-2 



U 
(72 

t' 



(All) 



such that the dressed locator Goo,<t {t " t') can be cast into 



tl ^J^-%t2 



G{)0,a{t — t') — • 



= O-C 
t' t 



=^0 

(7 t' 



O-C 

t 



1 



0-1 



3-0 

t' 



0-1 

O-C 

t 



t2 t3 
1^1 0-2 1 



0'2 



t4 

a-2 

3-0 
t' 



(A12) 



Although the cluster's undressed one-particle Green function 
O — I ^ I — O can be calculated exactly using the Lanczos 
method, the undressed cluster, n-particle Green function, 

12 n 

O — I ^ ^ ^ ^ ^ I — o ' is unknown for n > 1. 
It is clear that these functions cannot be calculated directly 
using the prescription provided by Wick's theorem because 
they include many-body Coulomb contributions coming from 
the impurity. 

In order to sort out this difficulty we propose another per- 
turbation expansion, assuming the cluster without the many- 
body term at the impurity as the unperturbed Hamiltonian and 
Hmb as the perturbation. The enormous advantage of this 
expansion in contrast with the previous one is that Wick's the- 
orem is applicable because the non-perturbed system is rep- 
resented by a one-particle Hamiltonian. For an infinite sys- 
tem, this expansion has been extensively used to calculate the 



one-particle Green function to study, for instance, the Kondo 
effect. In most cases, these studies have been restricted to ex- 
pansions in the self-energy up to second order in the Coulomb 
interaction parameter However, we are in a different situ- 
ation here because the system is finite and, more importantly, 
it requires the calculation of the Green function to all orders 
in the number of particles. In our case, the one particle Green 
function can be numerically calculated. After these diagrams 
are obtained, they are incorporated into the original diagram- 
matic expansion, Eq. ( IA12l i. in order to obtain the Green func- 
tion of the complete system Goo,o-(i — t'). When calculating 
the self-energy, this procedure in principle permits to sum up, 
to all orders in U, the most important families of diagrams. 
These are chosen among the ones that are essential to give a 
proper account of the region near the Fermi level. 

We use Eqs. (lAlb and (IA2l i to obtain this new diagrammatic 
expansion. It is worth mentioning that now the mean values 
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(• • • )o are calculated in the ground state of the cluster with- 
out the Coulomb interaction and that the evolution operator 
Eq. (IA2b requires the substitution of Hp by Hmb- 

In order to clarify the procedure and to establish the dia- 
grammatic rules, we calculate the first diagrams correspond- 
ing to the locator go,N,a,<yi (i, ^1,^2, t'), Eq. dATal l. We define 
three undressed Green functions. 



9oLit-t') = (T{co.(t),co.(0})o = O 



t' 

-o 



t 

o 



(A13a) 
a t' 



t 



(A13b) 
a t' 



(A13c) 

that, together with Eq. ( IA8I ), constitute the building blocks of 
the diagrammatic expansion. 

The contribution to the Green function to zero order in U, 
9o,N,a,cri {t, ti,t2: ^')' defined in Eq. (|A7at is given by 

50W,^(^:^l)30*W,<T(*2,t')'5<T<T' + 9ol]ait^'t')9NN,a'{tl,t2) = 

h , (A14) 

►t2 • ^ » 



tl 

a 

t 



a + 

6f' 



t 

o 



t' 
-o ■ 



From this result we infer that the contribution to Gqo cr(* 



t') in zero order in U is 



a 

t O 



t2 

a 

6t' 



h a' h 
a' 



The second diagram is a non-connected one and, as usual, 
does not contribute to the Green functioni^ 

To first order in U, incorporating all the possible contrac- 
tions resulting from the application of Wick's theorem and 
eliminating the non-connected diagrams, the contributions to 



t 



t' 



t 6 



t2 

a 

a 
Of. 



From these calculations, we conclude that there are two dif- 
ferent types of vertices • and o . At each • vertex there is 
one incoming and outgoing propagator and a factor of V has 



to be included. These are the vertices that result from the one 
particle Hamiltonian H^. The other vertex o comes from the 
Hamiltonian Hm b ■ There are two incoming and two outgo- 
ing spin CT and tr propagators and a factor of I] included at this 
vertex. As usual, the integral over the time variable associated 
to each vertex has to be taken. 

These rules are schematically represented as 



a 



To second order in U the topologically different connected 
diagrams that contribute to GoQ,cr{t — t') are 







(A15) 



The one-particle cluster Green function exactly obtained by 
numerical means, defined in Eq. dASb . can be thought of to be 
the result of the sum of the following infinite series of dia- 
grams: 




(A16) 



where (8) can be any site within the cluster although we are 
particularly interested in the impurity site o or the site • at 
the edge of the cluster. We use the dressed one-particle cluster 
Green function to incorporate all the diagrams of Eq. (IA16I ) 
into the expansion for the Green function 6*00,(7 ^ t'), 
Eq. iA5i . This results in 
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(A17) 



After taking a Fourier transformation in time, we define the 
self-energies, 

j:%{iu) ='i^^-^=gN+iA^)V^ (A18) 




where we have explicitly drawn the contribution to the self- 
energy up to terms proportional to U'^. 

The Green function of the system at the impurity can be 
written as a general Dyson equation: 

(A20) 

where i is restricted to be either or and the self-energy 
S°'(tj) is defined as 

muj) = ^%{uj)dm + So(w)'5.o. (A21) 

In order to compare the relative contribution of Sq and S^, 
which crucially depends upon the cluster size N, we proceed 
as follows. Considering that V ^ and using Eqs. ( IA18I ) 
and ( fT2] i we have that 

E^(c<j) oc 1/2 _ A-^^-i\ (A22) 

where we have ignored the Green function g7v+i,<T(w) be- 
cause we have numerically verified that in the neighborhood 
of the Fermi energy this function is independent of N. 

To evaluate the dependence of Eq upon N, we observe that 
all terms in Eq. ( |A19l l are multiplied by the square of the non- 
diagonal cluster propagator 

In addition, the dependence of this propagator on A is given 
by 

M{uj)UloU - A-(^-i)^/4, (A24) 
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where the function M{uj) goes asymptotically to zero when 
N increases above a characteristic length, which in our case 
corresponds to the size of the Kondo cloud. Defining f{N) ~ 
A'P{lu) we obtain, 

E^{uj) - f{N) X~iN-i)N/2^ (^25) 
As discussed in the main text in Sec. Illl the contribution to 



the self-energy Sq (w) can be neglected when compared with 
when the density of states of the leads is logarithmi- 
cally discretized, as their ratio is then proportional to: 



In this case 



^ f{N) X~i^~m2~i) . (A26) 
the embedding process is extremely simplified. 
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